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Abstract 

Observers at rest in a stationary spacetime flat at infinity can measure small amounts 
of rest-mass+internal energies+kinetic energies+pvessuve energy in a small volume of 
fluid attached to a local inertial frame. The sum of these small amounts is the total 
"matter energy" Em for those observers. If Mc 2 is the total mass-energy, the difference 
Mc 2 — Em is the binding gravitational energy. 

Misner, Thorne and Wheeler evaluated the gravitational energy of a spherically sym- 
metric static spacetime. Here we show how to calculate gravitational energy in any static 
and stationary spacetime with isolated sources with a set of observers at rest. 
, The result of MTW is recovered and we find that electromagnetic and gravitational 

3-covariant energy densities in conformastatic spacetimes are of opposite signs. Various 
examples suggest that gravitational energy is negative in spacetimes with special symme- 
tries or when the energy-momentum tensor satisfies usual energy conditions. 
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1 Introduction 



The paper deals mainly with stationary spacetimes of isolated sources that are asymptotically 
flat. 

In classical physics, energy has different forms which are additive. If energy is conserved 
one can assess how much change and transfers of energy occurs between the various forms. 
In Einstein's theory of gravitation different forms of matter energy are mixed in inseparable 
ways with gravitational binding energy. Gravitational energy is diffuse, partly mixed up 
with other forms, with matter energy, partly stored in the gravitational field itself. Alas, the 
matter is the source of gravity so that it is impossible to disentangle gravitational energy 
from other forms of energy with or without solving the field equations. This is a great loss 
that adds to that of "gravitational force" and "gravitational energy density" with which it 
shares a common origin: the principle of equivalence. 

One may, however, define gravitational energy with respect to a set of observers. Take for 
instance the matter-energy tensor of a perfect fluid; in standard notations 1 T£ = (pc 2 + 
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1 Indices A, fx, u, p, ■ ■ ■ = 0, 1, 2, 3; indices k,l,m,n ■ ■ ■ — 1,2, 3, the metric <? M „ has signature H and g 

is its determinant. Multiplication by y/—g is indicated by an additional hat like X on a previously defined 
symbol like X. Covariant derivatives are represented by a D and partial derivatives by a d. And for once 
neither G nor c are set equal to 1 and x° = ct. Finally the Levi-Civita symbol is e M „ pCT with £0123 = 1. 
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P)u IJ 'u u — 5^P and let be the 4- velocity field of local observers (g^ u w IM w L ' = 1)- If 
dV^ = j| e^ypudx 11 A dx p A dx a is the 3-dimensional coordinate volume element of a spacelike 
hypersurface [say, for instance, d 3 x on t = in adapted coordinates (x° = ct, x k )], the scalar 
element 

dE M = f»w v dV^ = + dE i > d(M c 2 )=pc 2 ^u^dV lx and u^" = -=L=. 

yl — v VI -f 

(1-1) 

d(Moc 2 ) is the element of rest mass + internal energy of the source, v the velocity of the 
matter with respect to the local observer in a local inertial frame and dEj is the external 
potential energy of the pressure on the small proper volume element ^/^gu^dV^. 

One "natural" selection of observers in static spacetimes which have a field of timelike 
Killing vectors of translations £^ are observers with velocities w^ 1 = They are at rest 

in any system of coordinates that singles out the inertial time at infinity £^ = {1,0,0,0}. 
They are also at rest with respect to the matter. There is, however, no such well defined set 
of observers in stationary spacetimes. For instance, in a Kerr spacetime one might prefer to 
choose the family of "zero angular momentum observers", the ZAMOS of Bardeen [3]. It is 
of course desirable to have some unique way of choosing observers in stationary spaces. Such 
a qualifying choice is not considered in this work. 

The sum Em of dEM's can be regarded as the total" rest mass-energy+internal ener- 
gies+kinetic energy + pressure energy" to which it reduces in the weak field limit. No grav- 
itational energy is involved in dEM- Therefore if Mc 2 is the total mass-energy of spacetime 
the difference Eq = Mc 2 — Em can be regarded as the total gravitational energy of spacetime 
for our set of observers. When gravity is a binding force, and this is not necessarily the case 
in general relativity, one expects that like in Newtonian gravity Eq is negative. 

Misner, Thorne and Wheeler [17] evaluated Eq along precisely this line of taught for 
spacetimes created by static spherical "stars" . Here we extend the calculation of Eq to any 
static spacetime with localized sources. The formula holds also in stationary spacetimes but 
its application may generate less enthusiasm because, as we said, sets of reference observers 
are not as well defined. In the few examples we treat Eq < when spacetimes have spe- 
cial symmetries or when the energy-momentum tensors of the source satisfy usual energy 
conditions. 



2 Conservation laws and observers 

(i) The conserved current for energy 

To begin with we need a conserved covariant vector involving the source of gravity Tjf. 
There are a number of available vectors, see in particular Abbott and Deser [1], Grishchuk, 
Petrov and Popova [7] or [10] 2 . All conserved currents for energy have the same form in 
stationary spacetimes: 

J M = (7? + (2-1) 

2 The Ashtekar and Hansen [2] conformal mass formula, the Hamiltonian formulation of Regge and Teitel- 
boim [21] and various formulas assembled by Szabados [24] to calculate energy in a finite volume (quasi-local 
energy) are integrals on closed surfaces. The integrant is an antisymmetric tensor whose diverge is a conserved 
vector which must be of the same form as eq. (2.1). One may thus start from those formulas as well but this 
needs a bit of additional work. 
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t^ u , which depends on the metric and its derivatives, is associated with part 3 of G-energy. 
#*„ differs from author to author reflecting the ambiguity in defining G-energy density. 

One necessary ingredient of J^ 1 is a background metric. In asymptotically flat spacetimes it 
is advantageous to introduce a flat background. The background is an extension from infinity 
inwards. It is an artifact to get a covariant description and is very useful when we want to 
use for instance spherical coordinates instead of Minkowski coordinates. Mappings on the 
background may be chosen at our convenience. The timelike Killing vectors associated with 
conserved energy are translations in the background, i.e. 

= (2.2) 

An overbar refers to covariant derivatives in the background. 

All covariant conserved vectors are also divergences of antisymmetric tensors or superpo- 
tentials J» v : 

> = Qjw (2.3) 

The superpotentials are not uniquely defined nor are they the same from author to author. 
However, the surface integral on a sphere at infinity of any valid superpotential must be the 
same i.e. Mc 2 . 

We shall here adopt 4 the current given in [9] and more fully described in [10]. There 
is a uniqueness [5] [8] to the "KBL-superpotential" which has been found satisfactory in all 
applications at spatial as well as at null infinity, on flat or non flat backgrounds including in 
the calculation of the very exotic mass of Kerr black holes on anti-de Sitter backgrounds in 
D-dimensions [6]. And it is quadratic in first order derivatives of the metric which may 
be a real boon as we shall see. In [10] t^ u is the covariant form of the quadratic Einstein 
pseudotensor. 

(ii) The KBL conserved current 

Formulations in [1], [7] and [10] are all on generally curved backgrounds. We shall start 
from the KBL identity I M = d v J> 11 ' which holds not only for curved backgrounds but also for 
arbitrary vector fields which are still denoted by £^ but are not necessarily Killing vectors. 
Our in (2.3) is a particular value of 1^. We shall first lay out the mathematical expression 
of and J^ v , which are somewhat complicated, and remind their physical content next. 
7 M = d v J^ u is originally a Noether identity which turns into a real conservation law once 
Einstein's equations G£ = kT^ are taken into account, n = ^f-. 

The KBL current and superpotential were derived in terms of Christoffel symbols of 
both the spacetime and the background Y or, more precisely, in terms of a tensor 
A^, = T^j, — We found it somewhat simpler to write J M = d v J^ v in terms of 

Lf = ^=D p €r , A, = l^L^ = ^lo gv ^ and F» = LJ» ' = D v gT (2.4) 
V 9 ^ 

A^„ and L '"'are so related: 

^ = " Vf + \ ( L V " 9 P a^) ■ (2.5) 

3 It is well known, see for instance [12], that the integral of Tj? dV^ in the weak field limit contains twice 
the Newtonian gravitational energy. 

4 The superpotential given in [1] does not provide the Bondi mass at null infinity [19] and the t^ v of [7] also 
given in [19] contains second order derivatives of the metric which is not a priori favorable to analyze the sign 
of the integrant. 
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The KBL superpotential which, in its original form is 

= Id^"] + If ^jfe"] with k u = -^=D^-gg^), (2.6) 



K 

can be rewritten as follows 5 : 



JV = - (P^D p ^ - L [H A e A + ^F u] >) where = - gW . (2.7) 
The current density I p is linear in f M , -D p £ M and covariant derivatives of 6 



= g + D °& in whi ch = g pX £ X (2. 



•ZpCT = if £ M is a Killing vector of the background. Thus 7 . 

J" = (f p -W)Z"+ (Jj p& R P *6!Z + K) C + ^\D~ v e + e upa D u z pW . (2.9) 
The three undefined quantities in this vector are 



t"„ = S»„-±6» V S with S% = ^ 



I L ppa - L pap \ L vpa - A P A U 



I1U 

* A 



(2.10) 

= _L (l x pv - F P 5 U X - 2L [H A ) , q pvpa = — {l pv g pa + l pa g pu - 2g^ p F>) . 

(2.11) 



The most significant properties of the conserved vector (2.9) were pointed out in [10]: (1) On a 
flat background (R X upa = 0) with a timelike Killing vector of translations, (2.2), equation (2.9) 
reduces to (2.1) where J p is the special value of I p in this case. The presence in the right 
hand side of T p ^ v shows that we are dealing with conservation of energy. The total energy 
is given by a surface integral of the superpotential (2.7). (2) If the Killing vector is one of 
rotations the surface integral gives the corresponding total angular momentum components. 
(3) On other backgrounds with various symmetries we obtain other conserved quantities. 

A remarkable feature of the conserved vector is that it is not useful for the very reason 
indicated in the introduction. There are no coordinate independent and background indepen- 
dent properties which one can associate with the different terms of J p . This will, however, 
change if we introduce a preferred set of local observers. 

(ii) Singling out the matter energy vector 

What we really want to see in J p is T p w u rather than T p £ v . We thus rewrite the first 
term of (2.9) like this: 

(f p -W) e = (f p -Wj w" + (f p -Wj r v where r v = - vf . (2.12) 



5 Our interest in the superpotential here is limited. For a detailed analysis of its structure and properties 
the reader is referred to the original paper. 

6 Z prT in [10] is here 2z p ^. Also indices in [10] are displaced with g never with g^ v . Here we do not stick 
to that convention. If an indice is displaced with the metric of the background g pl/ the indice has an overhead 
bar like in X-p — g^X" . Otherwise we write X p = g pv X v . 

7 In its original form, see [10], the last two terms in (2.9) are grouped slightly differently; in the present 
notations we would have written it like this a rip ^d [p ^ + <r M(pF) z pF + q^^D^z^ = a Aplf] d [p ^ + C" and C M 
is given by eq. (2.36) in [10]. 
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Ultimately w u will be a field of timelike unit vectors and £ u a timelike Killing vector field 
but for the moment both are unspecified to leave open the possibility to use the formulas 
in non stationary spacetimes and spacetimes with non-flat backgrounds. We then replace 
t v in (2.12) by its expression in terms of the Ricci tensor R pa using Einstein's 
equations in the foreground as well as in the background. Thus we set respectively 



T»t v = - 

K 



R pa = vJWRpa and t£t" = W^ a R pa . (2.13) 



A bar over the ro's indicates that the g's have been replaced by g's. Next we use the expression 
of R pa in terms of D-covariant derivatives and A pa as in [10] 8 : 

R pa = D x A x pa - D p A x x + A x pa Al v - A r j jX A x a + R pa . (2.14) 

We replace the (T — T)r's in (2.12) by (2.13) using (2.14), change the tz7(DA)-terms into 
\D(wA) — (Dw) A]-terms and obtain after straightforward but slightly tedious calculations 
the following expression for a modified conserved current J M and superpotential J^ v defined 
as follows: 

j^v = j»v_\_ (_l[H t \ + T [»pv]\ = I (Ip\p-d p ^ - L [H x w x + vMPA . (2.15) 

Notice that J^ v is similar to J^ u ; the £'s have been replaced by w's in the last two terms 
but not in the first one. X M is very similar to i 7 * with different coefficients: 

T" = (fj? -Jf)w v + W + + t p \D u w x ) + i- (Lf - S^F") D v i x + e upa D u z pW , 

(2.16) 

in this 

t^ x = - fs^F^ - L [H A ) and = ^- (i pa R pt7 w^ + l pp R p(7 r a - puX r x \ . (2.17) 

(Hi) Gravitational energy in stationary spacetimes 

The very complicated current J M in (2.16) with its superpotential T^ v in (2.15) is consid- 
erably simpler on a flat background, Vf- = 0, with a timelike Killing vector of translations, 
~z P a = and D v ^ x = 0; we denote it then by J^: 



= -{-L [iiV \w x + w [ ^F u ^ and J» = d u J^ = f»w v + + iT x D u w x ) . 

(2.18) 

This new conserved current has the following properties. Let # be a field of observers at 
rest, i.e. 

w fl = %- with — >^ + 0»(-) for r^oo, (2.19) 
£ r 

where r — > oo is the radius of the infinite sphere in the background whose metric in spherical 
coordinates is 

ds 2 = c 2 dt 2 - dr 2 - r 2 (d9 2 + sin 2 9dip 2 ) = c 2 dt 2 - dr 2 - r 2 duj 2 . (2.20) 

8 The expression (2.14) given in [10] and to which we refer, for the facility of the reader since it is on the 
web, was, however, published in 1940 by Rosen [20]. 
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If is replaced by the Killing vector £^ in J^ v it is equal to the KBL superpotential. 
Therefore the surface integral of at infinity is equal to Mc 2 because of (2.19): 

/ \ J^dS^ = -<f \ (-"& u \e + &F^) dS^ = Mc 2 , dS^ u = \ e^dxOAdx*. 

(2.21) 

But the matter-energy 

E M = I f£w u dV^. (2.22) 

J oo 

Therefore in a stationary spacetime the gravitational field energy with respect to our set of 
observers is given by an integral whose integrant is quadratic in first order derivatives of the 
field: 

E G = Mc 2 -E M = I (#>" + ^ A AV) dV (2.23) 

J oo ' 

This is an integral that gives the total gravitational energy in a stationary spacetime with 
localized sources. The integrant is quadratic in first order derivatives and is covariant. 

Now we examine the gravitational energy in some particular static and stationary space- 
times. 



3 Spherically symmetric static spacetimes 

Considered in isotropic coordinates, the metric is 

ds 2 = a 2 (dx ) 2 -b 2 dr 2 with dr 2 =Y,{dx k f. (3.1) 

k 

a and b are functions of r with a(oo) = b(oo) = 1. The background metric, is in Minkowski 
coordinates, 

ds 1 = (dx ) 2 - df 2 (3.2) 
The components of and the non-zero L's are as follows 

k3 k 

g 00 = — and g kl = -ab5 kl . Moreover with n k = — , (3.3) 
a r 

L fc oo = a 2 (3^ ~^jn k , L kmn = -b 2 + ^ n k 5 mn and (3.4) 

-fcOO _ 1 ( cfi a \^k T kmn _ 1 , a \ k x mn 



L fcUU = - ir---\n k , L« mn = -{- + - \n k b mn . (3.5) 
a z b^ \ b a J o 4 \b a J 

With these expressions one can calculate the gravitational energy: 

Eq = —— f (b-^dV where dV = b 3 d 3 x (3.6) 
K Joo v ' 

Here and later on, dV represents the proper volume element. This simple formula shows that 
Eg < 0. Moreover it gives a nice 3-coordinate independent gravitational energy density: 

«-=4( 6-1 )' 2 < 3 - 7 > 



Using the Schwarzschild solution for empty space in isotropic coordinates, we also find that 

If fc-'YW-i^ 

K Jr<TM V ' 2 r M 
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vm is the radial isotropic coordinate boundary of the matter. If the matter is a thin hollow 
shell the integral is zero. 

By definition Eq is the same as the ft of MTW [17]. To make the connection explicit 
consider the same metric in Schwarzschild coordinates 

ds 2 = a 2 dt 2 - B 2 dR 2 - R 2 du? where 6-^ = B and br = R. (3.9) 

dR 

Consider also one of Einstein's equations in isotropic coordinates Gq = kTq or, see [25]: 

6 3 r 2 

One easily find, using (3.9), (3.10) and the well known solution 9 of (3.10) for p = i.e. 



r 



(Vby}' = npc 2 . (3.10) 



+ (3.11) 
that (3.8) can be rewritten in the form given in [17] which in our notation reads: 

Eg = — j (l - -B" 1 ) P c 2 dV. (3.12) 
4 Eg for conformastationary metrics 

The gravitational energy can also be explicitly evaluated for the far less symmetric confor- 
mastationary metrics [22] which may be written as follows [11]: 

ds 2 = f(dx° - A k dx k ) 2 - r l dr 2 . (4.1) 

It has been noticed in [4] that this metric is in harmonic coordinates provided, here at least, 
that Ak is divergenceless. This leaves no gauge freedom except for global translations and 
Lorentz rotations. 

Solutions of Einstein-Maxwell's equations with this form of metric have been the object 
of much scrutiny in empty 10 and non-empty spaces [11]. The background metric is again 
c 2 dt 2 — dr 2 . The elements we need to calculate Eq are g^ u , and L^ v ' . First we write g^ 
and g» v : 

9oo = f , 901 = -fAi , 9m = -f~ 1 6 k i + fA k Ai , ^g~=f-\ (4.2) 
g 00 = r 2 -\A\ 2 , 5°' = -^ , g kl = -6 kl with \A\ 2 = Y,Al (4.3) 

k 

Next, following (2.19), 

™ = -t= , w k = 0. (4.4) 

We shall considerably simplify the formulas for the L's by using where possible a vector 
notation for 3- vectors with indices down 11 like A = {Ak}- Thus the non zero L's are 

L ooo = _fA.v(f- 2 -\A\ 2 ) , L k00 = -f 2 dk(f- 2 -\An (4.5) 
L ooi = fX.^Ai , L m = f 2 d k A u (4.6) 
Lkoo = ~2dkf , L m = f^dkipAi) , L kmn = -/ _1 9 fc (/ 2 ^l m Ai)- (4.7) 



which we shall only need for r — » co. 
10 See [22] for a revue of the subject. 
n with one exception see (4.14). 
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With these expressions one easily calculate E G which, as we shall see is negative definite. Set 
E G = [ e G dV with dV = (l - f 2 \A\ 2 ) 5 f~ld 3 x. (4.8) 
Then the gravitational energy density e G , a scalar in 3-space, is 

(1-/VT)" 1 • (4.9) 
To gain some insight into this formula let us consider a few special cases. 

(i) Static weak field 

This is not a particularly interesting case except to show that in the weak field approx- 
imation we recover the classical Newton G-energy. Indeed in this case ^4 = and, setting 
/ ~ 1 + 2cf)/c 2 , we have 

e G = -^IWI 2 . (4-10) 
4> is the Newtonian potential. This formula can also be recovered from (3.7) with (3.11). 

(ii) Static fields 

If A = but / is not nearly equal to one, the metric is that of Papapetrou [18] and 
Majumbar [16] and 

e G = ~^\£\ 2 where £ = -c 2 Vln^f. (4.11) 

£ acts as a potential force or what is also referred to as a "gravoelectric" force per unit 
mass on test particles (see below). Synge [23], who called those metrics "conformastatic" , 
showed that solutions of Einstein's equations of this form appear for static charged dust with 
electric and gravitational forces in exact balance; the electric charge density a = \[Gp. It is 
interesting to note that the electric field in this case 12 is E = V-\/7 an d the corresponding 
electromagnetic energy density 

e EM = ^~\E\ 2 = -e G . (4.12) 

87T 

There is no local force in the matter. There is also no local field energy density in this case. 
The total energy is equal to the mass-energy of the dust. All this seems to make good sense. 



£G = - 



Ak 



r'lV/l 2 + / 3 |V x A\ 2 + / 2 W ■ A x (V x A) 



(Hi) Non static fields 

In the general case, e G can be written in terms of the gravoelectric field £ and a "gravo- 
magnetic" field 

B = c 2 ^fV x A. (4.13) 

12 See for instance [15]. Incidentally on page 8 of that paper there is a printing mistake in the expression for 
E: there must be a 1/VC. 
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The reason for this electromagnetic analogy is that the equation for slow motions of test 
particles, in strong fields as seen in a 3-space 13 with metric f^ 1 5ki [14], with velocity v and 
momentum p, 



$={ v k = c <^ = ^} and p = fy f -l mv ky is tiL = m (e + -xB). (4.14) 
as dr ' dr \ c 



In terms of £ , B and A, expression (4.9) can be written 



1 



ec = - 



8vrG 



/(I - -J 2 \A\ 2 )\S\ 2 + \\B+ JfAx£\ 2 + ]f(A- £?] (1 - /VTr 1/2 - 



(4.15) 

One can see from this expression that ec is manifestly negative. For comparison we notice 
that electromagnetic energy density in this spacetime has a similar form. Denoting by E = 
{Ek = Fko = dkAo} and B = {Bk = £kmnd m A n } the electromagnetic vector fields: 



1 

8tt 



(1 - f 2 \A\ 2 )\E\ 2 + f 2 \B\ 2 + f 2 (A- E) 2 ] (1 - f 2 \A\ 2 Y 1/2 - 



(4.16) 



5 Eg of a weak stationary field 

In Minkowski coordinates in the background, the metric of the foreground is usually written 

ds 2 = {rjua, + h^dx^dx", (5.1) 
i]^ = diag (1, -1, -1, -1) and \h^ u \ < 1. Then, 

= _ _ 1 = _ ~ h »v with h = J]P a h ^ and L ^ ixv = _ dp lw ( 5 2 ) 

In harmonic coordinate, d v h^ v = 0, there is no further gauge freedom available except uniform 
translations and Lorentz rotations. Stationarity implies d^h^ = while Einstein's equations 
become 

V 2 h£ = 2kT{?. (5.3) 
Static observers have velocity components 

w° = 1 - - h° and w k = 0. (5.4) 

If we regard the hP k as components of a vector h in the background whose curl (V x h) = 
{J2 tmkidkh®}, and further decompose h™ into a traceless part t^C an d a trace h™, we find 
that 

eG = "i l ^ X ^ |2 "l^( l ^ |2 + ^ 1 ^ 12 " 2 ^™"^) where * = K-\K~hl 

(5.5) 

Note that the expression in terms of /i's instead of h's is exactly the same. We see that 
£g is not manifestly negative. However, if Vy/i™ = or is negligible, then G-energy is 
indeed negative. This is, for instance, the case if the source of gravity is a perfect fluid 
for which Tq = pc 2 and T™ = —P5™. In this case, Einstein's equations (5.3) imply that 

13 This is the 3-space in which proper lengths are measured. 
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V 2 (t^™) = because T™ is now traceless. From this follows, with usual boundary and 
regularity conditions, that Vrh™ = and, 

£G = __L|V x Hf - ^ (iVigP + i|V/,|| 2 ) < 0. (5.6) 

The equation of slow motion of a test particle in the weak field has the form (4.14) with 

_h0 _ _ 

£ = -c 2 V^ and B = c 2 V x h. (5.7) 

If |V^| <C |V/io|, like in a perfect fluid with non relativistic pressure or tension, \P\ <C pc 2 , 
the gravitational energy density is: 

^ = -^^T + 7l^l 2 V M) 



8ttG V 4' 

The source of Eq may be a fluid in fast, but not too fast and not necessarily uniform rotational 
motion, with vorticity, for instance, like in a Dedekind ellipsoid. 

Notice that Landau and Lifshitz [14] prefer to regard the gravomagnetic force as a Coriolis 
force produced in an orthogonal frame rotating with absolute angular velocity Q = |V x h. 



6 Comments 

In a similar vein we can calculate gravitational energy of cosmological perturbations on a 
Robertson- Walker background. These backgrounds all have a conformal timelike Killing field. 
Thus one may be facing "conformal matter-energy" and "conformal G-energy" . Conservation 
laws and additivity of different forms of energy remains valuable knowledge but the conserved 
vector may not have the simple form it has here in (2.17) and the interpretation of globally 
conserved quantities is quite different from asymptotically flat spacetimes. 

It is possible that Eq is negative for any stationary metric with reasonable physical 
sources, that is, those whose energy- momentum tensor satisfy reasonable energy conditions. 
It is, however, interesting to notice that Eq is negative in some spacetimes with special 
symmetries irrespective of any energy condition. 
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